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source, ledit proc6d6 mettant en uvre une transformation en ondelettes associant o una image source une image sous-^chantillonnr 
^SZ'™ J::rd^' f^^^^^ ondelettes conespondan. 6 au moins une image de ditails. pour au 7Z „ nW^^Se 

position, un point d int&6t 6tant un point associ6 6 une region de l image pr&entant des hautes fi^uences Le proc6d6 comorend 1« 
stapes suivantes :- application de ladite transformation en ondelettes 6 ladite image source ;- constSrdi^e'^^s^Siie 
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Proc^d^ et dispositif de detection de points d'int^rgt dans une image 
num^rique source, programme d'ordinateur et support de donn^es 
correspondants. 

5 1. domaine de rinvention 

Le domaine de rinvention est celui de la detection de points d'interdt, 
encore appel6s points saillants dans vine image numerique. Plus pr6cis6ment, 
rinvention concerne une technique de detection de points d'int6r6t mettant en 
o&uvre une approche de ondelettes. 
10 Un point d'int6rgt pent gtre consid€r6 comme le iepr6sentant d'une region 

spatiale de I'image vehiculant une partie importante de Tinfomiation. 

Historiquement, la notion de point saillant a 6t6 propos6e dans le domaine 
de la Aosion par ordinateur, oft Tun des probl^mes majeurs consistait k d6tecter les 
coins des objets (d'oili le t^me « saillant », utilise par la suite k titre de synonyme 
15 du terme « d'int6r6t »). Plus tard, cette notion a 6t6 elargie k d'autres 
caracteristiques des images comme les contours, les jonctions, etc. 

Dans le domaine du traitement d'images, la detection des points saillants 
correspondant aux coins des objets n'a que pen d'int€r8t En effet, les coins sont 
g6n6ralement des points isol6s, ne repr6sentant qu*une faible partie de 
20 information contenue dans Timage. De plus, leur detection ginhrc des amas de 
points saillants dans le cas de regions textures ou bruit6es. 

DifKrentes autres techniques ont 6t6 proposfes, s'int^ressant notamment 
aux points saillants correspondant aux zones de hautes frequences, c'est-k-dire 
aux contours des objets. L'invention s'applique plus pr6cis6ment k ce type de 
25 technique. 

On pr^sente ci-aprfes plus en detail diff^rentes techniques connues de 
detection de points saillants. 
2, JEtatdeF^rt 

La detection de points saillants (aussi appel6s points d'int6r6t) dans les 
30 images est un problfeme qui a suscit^ de nombreuses recherches depuis plusieurs 
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anuses. On pr^sente dans cette section les principales approches retenues 
classiquement dans la Utt^rature. On pourra se r6f6rer m document [5] Qes 
documents auxquels fl est fait r^fdrence sont regroup6s en Annexe B) pour un 6m 
de Tart plus d6taiI16. 

5 Une des piemi&ies m^fliodes a 6t6 prpposfe par Harris et Stephens [7] pour 

la detection des coins. Ce type de points etaient alors consid6r^s comme 
vdhiculant une quantity importante de I'information et trouvait son appUcation 
dans le domaine de la vision par ordinateur. 

Pour d6finir ce d^tecteur, on calcule en chaque point p(x,y) de I'image /. la 
10 quantity : 



25 



oil est une matrice definie par : 



Il(x.y) IJx.y)I/x,y) 
/x(x.y)I/x.y) Il(x,y) 

oil : 

*♦* G( o) denote un noyau gaussien de variance ; 

❖ ® denote le produit de convolution ; 

❖ /, (resp. /,) denote la d&ivfe premiere de / suivant la direction x (resp. 

❖ £>cifM^y) denote leddtetminantde la matrice M • 
20 ❖ ^Af*,) denote la trace de la matrice Af^; 

❖ k est une constante g6n6alement utilise avec une valeur de 0.04. 

Les points saiUants sont alors d6finis par les extrema locaux positif de la 
quantity 

Dans [5], les auteurs proposent aussi une version plus precise du d^tecteur 
de Harris et Stephens. Cette version remplace le calcul des d&iv^ de I'image / 
par un calcul precis des d^rivdes du noyau gaussien. 

Le d^tecteur de Harris et Stephens pr&ent6 ci-dessus a 6X6 ^tendu au cas 
des images couleurs dans [6]. Pour ce faire. les auteurs ^tendent la ddfinition de la 
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matrice qui devient alors : 



10 



15 



20 



oil : 



❖ R^G^B^ d6notent respectivement les d6rivdes premieres des plans 
colorimdtriques rouge, yett et bleu dans la direction x ; 

❖ Ry.Gy,By d^notent respectivement les ddriv^s premieres des plans 
colorimdtriques rouge, vert et bleu dans la direction y ; 

Dans [10], les auteurs considferent les points saillants conune les points de 
I'image pr&entant un fort contraste. Pour constniire un tel d^tecteur, les auteurs 
utilisent une approche multi-resolution bas6e sur la construction d'une pyramide 
gaussienne. 

Supposons que I'image / soit de taille IV. Nous pouvons d^finir une 
pyramide a N niveaux oil le niveau 0 correspond k I'image originale et le niveau 
N-1 correspond ^ une image de 1 pixel. 

Au niveau k de la pyramide, le contraste au point P est d^fini par : 
^k(P) °- ^*^pj avecOg ksN-let C^(P) = 1 

oil G^P) d6finit la luminance locale au point P et au niveau k, et B^P) d^finit la 
luminance du fond local au point P et au niveau k. 

Ces deux grandeurs sont calcul6es en chaque point et pour chaque 
niveau de la pyramide. On pent done les repi^ senter par deux pyramides appel6es 
pyramide de luminance et pyramide du fond et d^finies par : 

y,'^(M)G^_/M) 

B,(P)^ 2W(Q)G,,,(Q) 

QeParait(P) 

oh. : 

❖ Les notations Fils(P) et Parent(P) ddnotent les relations de hi&archie 
25 dans la pyramide gaussienne ; 



wo 2004/068410 



4 



PCT/FR2003/000834 



❖ w est une fonctibn de poids nonnaUsde pouvant fitre i^gl^e afin de 
simiileT la pyramide gaussienne ; 

❖ W est une fonction de poids normalis6e prenant en compte la fagon 
dont P est utilise pour constniue la luminance de ses ascendants dans 

5 la pyramide. 

Dans cette approche, un point saillant est un point caract6ris6 par une 
valeur 61ev€e du contraste local. Afin de prendre en compte la non sym6trie de la 
grandeur Q, les auteurs introduisent une nouvelle grandeur afin d'obtenir une 
valeur nuUe pour une situation de non contraste et une valeur sup6rieure k z6ro 
10 partout ailleurs. 

Cette nouvelle grandeur est d^finie par : 

cKP) = MJ^J^<lhMEA WltiMlA\ 

\ B,(P) ' 255 -B,(P) ]' 

Avec cette nouveUe grandeur, les points saillants sont d^finis par les 
maxima locaux de Cl sup6rieurs ^ un seuil fix6. 

Le detecteur de points saillants pr&entfi initialement dans [11] est sans 
dome le plus proche de la pr^sente invention puisqu'il est aussi bas6 sur 
I'utilisation de la throne des ondelettes. En effet. les auteurs considferent que les 
points v6hiculant une partie importante de rinformation sont locaUsds dans les 
regions de I'image prSsentant de hautes frequences. 

En utilisant des ondelettes k support compact, les auteurs sont capables de 
determiner un ensemble de points du signal /(que I'on suppose pour I'instant 
mono-dimensionnel) qui ont 6te utiUs^s pour calculer n'importe quel coefficient 
ondelette D^jf( n) , et ce, ^ n'importe queUe resolution 2^ (j & -1) . 

A partir de cette constatation. une hierarchie de coefficients ondelettes est 
construite. Cette hierarchic determine pour chaque niveau de resolution 2^etpour 
chaque coefficient ondelette D^, f(„)d<i ce niveau. I'ensemble des coefficients 

ondelettes du niveau de resolution immediatement sup&leur 2^*'necessaiies pour 
calculer D^,f(n) : *^ 
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CiD^tfin)) = [d^^ f(k),2n &k£2n + 2p-l},0sn<2JN 

ot p denote la r6gulaiit6 de la base d'ondelette utilis6e (i.e. la taille du filtre 
ondelette) et N d6aato la longueur du signal original /. 

Ainsi. chaque coefficient ondelette D^jf(n)ost calcul6 k partir de 2-Jp 
points du signal/ Ses coefficients fils C(D^J(n)) donnent la variation d'un 
sous-ensemble de ces 2'^ p points. Le sous-ensemble le plus saillant est celui 
dont le coefficient ondelette est maximal (en valeur absolue) au niveau de 
resolution 2-'*'. 

II convient done de consid&er ce coefficient k ce niveau de resolution. En 
apphquant de fa9on recursive ce processus, un coefficient D^^f(n)est 

s61ectionn6 k la resolution 1 Ce coefficient repr^sente 2p points du signal/ Pour 

seiectionner le point saillant correspondant dans/ les auteurs proposent de retenir 
parmi ces 2p points, celui dont le gradient est maximal en valeur absolue. 

Pour etendre cette approche aux signaux bi-dimensionnels que sont les 
images, tes auteurs ^pUquent la mdme approche k chacune des ttois sous-bandes 
D^jI.D^jI,D^jQ^ I denote I'image originale. Dans le cas des images, le support 
spatial de la base d'ondelette est de taille 2pi,2p. Ainsi. le cardinal de 
C( ly^J (x,y)) est de 4f pour tout *=7,2,5. Pour chaque orientation (horizontale, 
verticale et oblique), la methode cherche. parmi les coefficients ffls d'un 
coefficient domie. celui dont I'amplitude est maximale. Si differents coefficients 
de differentes orientations amfenent au mSme pixel de /, alors ce pixel est 
consid^ comme un point saillant 

Cette technique a ete notamment utilise en indexation d'images dans [9]. 

Inconvenients des teelii»<n ^es ant^riMirP^ 

Comme cela est montre dans la section precedente, de nombreuses 
meihodes ont ete proposees dans la Utteratuie pour la detection de points saillants. 

La di£f(&:ence majeure entre ces approches repose sur la definition m6me 
de point saiUant Historiquement. les chercheurs dans le domaine de la vision par 
ordinateur s'interessaient aux coins des objets. C'est ainsi qu'a ete propose le 
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detecteur de Harris et Stephens [7] qui a 6t6 r^nunent 6tendu k la couleur dans 
[6]. Les coins des objets ne leprfeentent cependant pas une information pertinente 
dans le domaine du Iraitement d'image. En effet, dans le cas d'images faiblement 
textur6es, ces points vont Stre dpaipill6s dans I'espace et ne donneront pas une 
5 representation satisfaisante de rimage. Dans le cas d'images textures ou bruitfes. 
les points saillants seront tons concentres dans les textures et donneront une 
representation locale et non globale de I'image. 

La definition de la saiUance basde sur le contraste [10] est nettement plus 
interessante pour le traitement d'image. Malheureusement, cette approche souffie 
10 du meme d^faut que la pr^c^dente dans le cas de regions textur&s ou bruits. 

❖ L'approche basde sur les ondelettes proposde par E. Loupias et N. Sebe 
[11] est manifestement la plus robuste et la plus interessante. En effet, on 
salt depuis longtemps que les contours repr^sentent I'information 
primordiale d'une image puisqu'eDe est en parfaite adequation avec le 
15 systfeme visuel humain. 

^ Obiectifs et caract^rifiH gucs l»inY «'"«»" 

L'invention a done notamment pour objectif de paUier les differents 
inconvenients de I'etat de Tart. 

Plus precisement, un objectif de l'invention est de foumir une technique de 
detection de points saillants conespondant k une haute frequence, et ne 
privaegiant aucune direction particuKfere dans I'image. 

Un autre objectif de l'invention est de foumir une telle technique, qui 
necessite un nombre r^ d'opetations. par rapport aux techniques comiues. 

Notamment, un objectif de l'invention est de foumir une teUe technique 
permettant d'utiliser des bases d'ondelettes avec un support de taiUe importante. 

Ces objectifs, ainsi que d'autres qui apparaltront plus clairement par la 
suite, sont atteints k I'aide d'un procede de detection de pomts d'int^t dans une 
image nmnerique source, ledit procede mettant en oeuvie une transformation en 
ondelettes associant k mie image source une image sous-^chantiUonnee. dite 
image echeUe. et des coefficients ondelettes conespondant k au moms une image 
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de details, pour an moins un niveau de decomposition, un point d'int&6t 6tant un 
point associd k une region de I'image pi^sentant des hautes ft^uences. 
Selon I'invention, ce proc6d6 comprend les Stapes suivantes : 

- appUcation de ladite transformation en ondelettes h ladite image 
source ; 

- constriction d'une arborescence unique, k partir des coefficients 
ondelettes de chacune desdites images de details ; 

- selection d'au moins un point d'intdrSt, par analyse de ladite 
arborescence. 

Dans le present document, on appelle par simplification image source une 
image d'origine ou une image ayant subi un pr6traitement (calcul de gradient, 
changement d'espace colorim^trique,. . .). 

De fagon avantageuse, on determine, pour chaque niveau de 
decomposition, au moins deux images de details coirespondant respectivement k 
15 au moins deux directions pr6d6tennin6es par ladite transformation en ondelettes. 

Cette transformation en ondelettes peut notamment utiliser des ondelettes 
de premiere ou de deuxifeme g6n6ration (basde maiUage). 

Notamment. lesdites images de d^taUs peuvent compiendre ; 

- une image de details repr^sentant les hautes ft^uences verticales ; 

- une image de details repr&entant les hautes ft^quences horizontales ; 

- «ne image de details repr&entant les hautes ft^iuencesdiagonales, 
De fagon avantageuse, le proc6d6 de I'invention comprend une 6tape de 

fusion des coefficients desdites images de details, de fagon k ne piivil6gier aucune 
direction de ladite image souice. 

De fagon avantageuse, ladite 6tape de consti^ction d'une arborescence 
repose sur une approche de type Zerotree. 

Ainsi, pr6f€rentiellement, chaque point de I'image dcheUe de resolution 
minimum est la racine d'un arbre k laqueUe on assode un noeud fils k chacun des 
coefficients ondelettes dechacune de la ou desdites images de details localises au 
m6me emplacement, puis on associe r^cursivement, k chaque noeud fils d'un 
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niveau de resolution donn^, quatre noeuds fils form6s par les coefficients 
ondelettes de I'image de details de mSme type et du niveau de r&olution 
precedent, et associ^e h la region corcespondante de I'ima^ source. 

Selon un aspect avantageux de I'invention, ladite 6tape de selection met en 
5 ceuvre une 6tape de construction d*au moins une carte de saillance, affectant 
auxdits coefficients ondelettes une valeur de saillance representative de son 
int^rgt Prefgrentiellement, on construit une carte de saillance pour chacun desdits 
niveaux de r6solution. 

De fagon avantageuse, pour chacune desdites cartes de saillance, on 
10 fusionne pour chaque valeur de saiUance les informations associ6es aux trois 
coefficients ondelettes correspondant aux trois images de detail, de fagon k ne 
privil6gier aucune direction dans Timage. 

Selon un aspect prgf&entiel de I'invention, une valeur de saillance d'un 
coefficient ondelette donn6 d'un niveau de resolution donne prend en compte la 
15 ou les valeurs de saiUance des coefficients ondelettes descendant dans ladite 
arborescence dudit coefficient ondelette donn6. 

PreferentieUement. une valeur de saillance est une relation lineaire des 
coefficients ondelettes associ^s. 

Dans un mode de reaUsation particuUer de I'invention, la valeur de 
20 saillance d'un coefficient ondelette donne est calcuiee k partir des Equations 
suivantes : 

Dans ces equations, le param^tre pent par exemple valoir -1/r pour 
toutes les valeurs de k. 

25 Selon un autre aspect pteferentiel de I'invention, ladite etape de selection 

comprend une etape de construction d'une arborescence desdites valeurs de 
saillance, repose avantageusement sur une approche de type Zerotree. 
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Dans ce cas, ladite 6tape de sflection comprend avantageusement les 
Stapes de : 

- tri par ardre d^croissant des valeurs de saiUance de la carte de saiUance 
correspondant k la resolution miTiitniiTn ; 

5 - selection de la hranchepr6sentant la valevirde saiUance la plus eiev^e pour 
chacun des artnes ainsi tci6s. 

Selon un aspect pr^f&entiel de I'invention, ladite etape de s^ection de la 
branche pr^sentant la valeur de saiUance la plus 61ev6e met en oeuvre un parcours 
de I'arbre correspondant k partir de sa racine, et une selection k chaque niveau de 
10 I'arbre du noeud ffls pr^sentant la valeur de saiUance la plus 6lev6e. 

Comme d^j^ mentionn6, I'invention permet d'utUiser de nombreuses 
transformations en ondelettes. Dans un mode de realisation particuUer, on met en 
oeuvre la base de Haar. 

Dans un mode de r6aUsation particuUer, on choisit comme niveau minimal 
15 de resolution 2"*. 

Le precede de I'invention peut par aiUeurs comprendre une etape de calcul 
d'une signature d'image. k partir d'un nombre piedetermine de points d'int6rgt de 
ladite image. 

Ladite signature peut ainsi notamment Btxe utilisee pour I'indexation 
20 d'images par leur contenu. 

Plus generalement I'invention tix)uve des appUcations dans de nombreux 
domaines, et par exemple pour : 

- le tatouage d'images ; 

- I'indexation d'images ; 

25 - la detection de visage(s) dans une image. 

L'invention conceme egalement les dispositifs de detection de points 
d'interet dans une image numerique source mettant en oeuvre le precede tel que 
decrit ci-dessus. 

L'invention conceme encore les programmes d'ordinateur comprenant des 
instructions de code de programme pour I'execution des etapes du procede de 
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detection de points d'int^rgt d6crit ci-dessus, et les supports de donn6es 
num&iques utilisables par un oidinateur portant un tel programme. 

D'autres caract^stiques et avantages de I'invention apparaitront k la 
lecture de la description suivante d*un mode de realisation pr^fdrentiel, donn6 h 
5 title de simple exemple illustratif et non limitatif, et des dessins annexes parmi 
lesquels : 

la figure 1 illustre le principe de I'analyse multi-rdsolution d'une 
image I par transformation en ondelettes ; 

la figure 2 pr6sente de fa9on sch6matique une transformation 
10 ondelette ; 

la figure 3 est une representation d'une arborescence de 
coefficients ondelettes selon I'invention ; 

la figure 4 prfisenle un exemple de cartes de saillance. et les arhies 
de saillance correspondants ; 

la figure 5 illustre la saillance d'une branche de I'aibre de la figure 
4; 

les figures 6a et 6b iUustrent des resultats experimentaux du 
procede de I'invention. la figure 6a pr6sentant deux images 
originales et la figure 6b les points saillants correspondants ; 
la figure 7 iUustre un precede d'indexation d'images mettant en 
oeuvre le precede de detection de I'invention. 
5. Identifleatinn des a^me nts technianes essenticls de Pinveiitinn 

5.0 Piinctnes g^nA^nwnr 
Un but de I'invention est done la detection des points saiUants d'une image 
/. Ces points correspondent aux pixels de / appartenant k des regions de haute 
frequence. Pour ce faire, on se base sur la theorie des ondelettes [1][2][3]. 
L'Annexe A presente une br&ve presentation de cette dieorie. 

La transfoimee en ondelettes est une representation multi-resolution de 

I'image qui permet d'exprimer I'image aux differentes resolutions - - etc 

2'4' 

30 Ainsi, k chaque niveau de resolution 2' (J s -I) . la transformee en ondelettes 



20 



25 



wo 2004/068410 



11 



PCT/FR2003/000834 



lepr&ente I'image /, de taille nxm'-l'' xl^ (k,l^Z) , sous la forme : 

❖ d'une image grossi^ie A^J ; 

❖ d'une image D\,I6& details iepi€sentant les hautes fc^uences verticales 

(ie. les contours horizontaux) ; 
5 ❖ d'une image Dlld& details repr€sentant les hautes frequences 

horizontales (i.e. les contours verticaux) ; 

❖ d'une image D^^/de details lepr^s^tant les hautes frequences diagonales 

(i.e. les coins). 

Chacune de ces images est de taille 2**-' x 2'*^ . La figure 1 illustre ce type 
10 de representation. 

Chacune de ces trois images est obtenue k partir A^j^,I par un filtrage suivi 

d'un sous-echantillonnage d'un facteur deux dans chaque direction comme le 
montre la 2. n est k noter que Ton a A^I = I. 

L'invention consiiste h choisir, en premier lieu, une base d'ondelettes et un 
15 niveau de resolution minimal 2'(r s -1) . Une fois la transformation ondelette 
effectuee, nous proposons de parcourir chacune des trois images de details D\l , 
£>J/ et Dj/afin de construire une arborescence de coefficients ondelettes. Cette 
arborescence est basee sur I'approche Zerotree [4], initialement proposee pour le 
codage d'image. Elle permet de mettre en place une carte de saillance de taille 

20 2**' x2'*' refletant I'importance de chaque coefficient ondelette k la resolution 
2Yrs-i;. 

Ainsi, im coefficient ayant une saillance importante correspond k une 
region de /presentant des hautes frequences. En effet. un coefficient ondelette de 
module important k la resolution 2Yr s -i; correspond k un contour de I'image 
25 Ay^I suivant une direction particuUdre (horizontale, verticale ou oblique). 

L'approche Zerotree nous indique que chacun des coefficients ondelettes k la 
resolution 2' correspond k une zone spatiale de taille 2"'' x2-''dans I'image /. 

A partir de la carte de saillance construite, l'invention propose une 
methode permettant de choisir parmi les 2"' xl"' pixels de /, le pixel le plus 
30 representatif de cette zone. 
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En termes d'applications potentielles, la detection de points saillants dans 
les images pent Stie uti]is6e, de fagon non exhaustive : 

❖ Pour le tatouage d'images. Dans ce cas, les points saillants donnent des 
indication quant k la localisation possible de la marque afin de garantir sa 
lobustesse ; 

❖ Pour Tindexation d'image. En d^tectant un nombre fixe de points saillants, 
on pent en d6duire une signature de I'image (bas6e, par exemple, sur la 
colorim6trie autour des points saillants) qui peut ensuite 6tre utilisde pour 
le calcul de similarity inter-images ; 

❖ Pour la detection de visages. Parmi les points saillants conespondant aux 
hautes fi:i6quences de Timage, certains d'entre eux sont localises sur les 
carac^ristiques faciales (yeux, nez, bouche) des visages presents dans 
I'image. Ds peuvent alors Stre utilises dans un ptoc6d6 de detection de 
visages dans les images. 

La technique de I'inv^tion difffere de celle propos^e par E. Loupias et N. 
Sebe [1 1]. Les principales differences sont les suivantes : 

❖ L'algorithme de recherche des points saillants propos6 par Loupias et Sebe 
requiert une recherche parmi 2^^' x 4p^ x Scoefficients pour chaque niveau 
de resolution 2^et pour une image carr€e. Notre algorithme est 
ind^pendant de la taille du support de la base d'ondelette, ce qui nous 
amfene k une recherche parmi 2^-'x4x3coefficients. Get avantage nous 
pennet de pouvoir utiliser des bases d'ondelette avec un support pouvant 
6tre de taille inq>ortante alors que la plupart des publications utilisant le 
ddtecteur de Loupias et Sebe utiKsent la base de Haar, loin d'6tre optimale. 

La mdthode de Loupias et Sebe considerent les sous-bandes 
independamment les unes des autres ce qui les amfenent k detecter en 
priority les points de gradient maximal dans toutes les directions (i.e. les 
coins). De notre c6te, nous fusionnons Tinformation contenue dans les 
difKrentes sous-bandes ce qui nous permet de ne privil^gier aucune 
direction particuli^re. 
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S.1 Tranformation ondelette 
La transformation ondelette est un outil math^matique puissant permettant 
I'analyse multi-r6solution d'une fonction [1I2][3]. L' Annexe A pour un rapide 
survol de cet outil. 

5 Dans le cas de Tinvention, les fonctions consid6i^es sont des images 

num6riques c'est-^-dire des fonctions bi-dimensionnelles discretes. Sans perte de 
g6n6ralites, nous supposons que les images trait^s sont 6chantiIlonnees sur une 
grille discrete de n lignes et m colonnes et k valeur dans im espace de luminance 
^chantillonnd k 256 valeurs. De plus, nous supposons que « = 2* (k&Z) et que 
10 m = 2'(lEiZ). 

Si nous notons / I'lmage originale, nous avons alors : 
^.|[0,w]x[0,»] — [0,255] 

Comme mentionnd k la section4, la transformation ondelette de / permet 
une representation multi-resolution de /. A chaque niveau de resolution 
15 2^(j^-l), la representation de / est donnee par une image grossidre A^j I et par 
trois images de details D\jl , 0^1 &t Dl,I. Chacune de ces images est de taille 

2**^ X 2'*^ . Ce processus est illustre sur la figure 2. 

La transformation ondelette necessite le choix d'une fonctiion echelle 
<b( x) ainsi que le choix d'une fonction ondelette ^(x). A partir de ces deux 
20 fonctions, on derive un filtre echelle H et un filtre ondelette G dont les reponses 
impulsionnelles lespectives h&tg sont definies par : 

hin) = {^^.. - n))Vn e Z 

8{n) = {v;^., («),^(m - n))V» S Z. 

Notons respectivement iret Gles filtres miroirs de H et G (i.e. 
h(n) = h(-n) et g(n) = g(-n)). 

25 On pent alors montter[l](cf. figure 2) que: 

^* ^ij/^P®^' calcuiee en convoluant 4jy./avec ^dans les deux 

dimensions et en sous-echantillonnant d'un facteur deux dans les deux 
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dimensions ; 

❖ D^y/peutetrecalculfeen: 

1. convoluant 4py*,/avec ffsuivant la direction y et en sous- 
6chantillonnant d'un facteur deux suivant cette meme direction ; 
5 2. convoluant le resultat de Fetape 1) avec G suivant la direction x et 

en sous-echantiUonnant d'un facteur deux suivant cette mSme 
direction. 

❖ D^y/peutStrecalculeeen: 

1. convoluant ^^y+i/avec G suivant la direction y tt en sous- 
10 &:hantillonnant d'un facteur deux suivant cette mSme direction ; 

2. convoluant le resultat de T^tape 1) avec H suivant la direction x et 
en sous-6chantillonnant d*un facteur deux suivant cette mSme 
direction. 

❖ iJ^^/peutStrecalculeeen: 

15 1- convoluant ^^^x/avec G suivant la direction y et en sous- 

6chantillonnant d'un facteur deux suivant cette mgme direction ; 
2. convoluant le resultat de I'etape 1) avec G suivant la direction x et 
en sous-6chantillonnant d'un facteur deux suivant cette mSme 
direction. 

^0 Construction de Vixrh orescence des coefficients ondelettex 

Une fois la transformation ondelette effectu6e jusqu'^ la resolution 
!''( r !S -1^ , on dispose : 

❖ D'une image approximfe A^J ; 

❖ De trois images de details Dj,/, Dj/, 0^1 par niveau de resolution 
25 2^avecj=-l,...^. 

On construit alors une arborescence des coefficients ondelettes en se 
basant sur la technique Zerotree [4]. Les arbres sont construits de la fafon 
suivante (cffigure 3) : 

❖ Chaque pixol p(x,y) de I'image A^I est la racine d'un arbre ; 
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❖ Chaque racine p(x,y) se voit attribuer trois noeuds fils d6sign6s par les 
coefficients ondelettes des trois images de details D^/(s=lA3) localises 

au meme emplacement (x,y) ; 

❖ Du fait du sous-6chantillonnage d'un facteur deux effectu6 par la 
5 transformation ondelette & chaque changement de resolution, chaque 

coefficient ondelette a'y(x,y) (s=l,2,3) correspond k une zone de taille 

2x2pixels dans I'image detail correspondante a la resolution 2"*^. Cette 
zone est localis6e en (lx2y) et tons les coefficients ondelettes y 
appartenant deviennent les ncEuds fils de a*, (x,y ) . 

10 De fagon recursive, on construit Tarborescence dans laquelle chaque 

coefficient ondelette al,(x.y)is=\,2,3 et 0>u>r) poss^de quatre nceuds ffls 
d6signes par les coefficients ondelettes de I'image D'^I localises dans la i^gion 
situeeen (2x.2y) etde taille 2x2pixels. 

Une fois I'arborescence construite, chaque coefficient ondelette 
15 a* ( x.y) (s=l,2, 3) correspond k une r6gion de taille 2"' x2"''pixels dans I'image 
detail Dr,I. 

5.3 Cnnstr uction des cartes de saiUance 

A partir de I'arborescence obtenue par I'etape precedente, nous proposons 
de construire en ensemble de -r cartes de saillance (i.e. une carte de saillance par 
20 niveau de resolution). Chaque carte de saiUance S^Q=-l„..;[) reflate 

I'importance des coefficients ondelettes presents ^ la resolution correspondante 
Ainsi, plus un coefficient ondelette sera juge important au sens de 
rinformation qu'il vehicule, plus sa valeur de saillance sera importante. 

n est k noter que chaque coefficient ondelette priviiegie une direction 

25 (horizontale, verticale ou oblique) suivant I'image de details h. laquelle il 
appartient. Cependant, nous avons choisi de ne priviiegier aucune direction 
particuliere et nous avons done fiisionne I'information contenue dans les trois 
coefficients ondelettes alj(x.y).al(x,y).al(x,y)qvi&lque soit le niveau de 
resolution 2^ et quelque soit la localisation (x,y) avec 0 s x < 2**-' etOsy< 2'*^ . 

30 Chaque carte de saillance S^j est de taille 2**-' x 2'*^ . . 
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D'autre part, la saillance de chaque coefficient h la resolution 2^ doit 
prendre en compte la saillance de ses descendants dans Tarborescence des 
coefficients. 

Afin de prendre en compte Tensemble de ces propri6t6s, la saillance d'un 
5 coefficient localise en (x,y) k la resolution 2^ est donn^e par la relation recursive 
suivante : 



Equation 1: expression de la saillance d*un coefficient 



10 ou 



Max( D^j ) (s=l,2,3) denote la valeur maximale des coefficients ondelettes 
dans Titnage detail D'^jl ; 

❖ (O £ £; i;)permet de regler r importance des coefficients de saillance 
suivant le niveau de resolution. II est k noter que Ton a = 7 . 

15 ❖ n est important de noter que les valeurs de saillance sont normalisees i.e. 

Comme on pent le voir sur I'l^uation 1, la saillance d*un coefficient est 
une relation lineaire des coefficients ondelettes. En eflfet, comme nous I'avons 
mentionne en section4, nous considerons les points saillants comme des pixels de 
20 rimage appartenant k des regions de haute frequence. Or, un coefficient ondelette 
<^l}(x»y)i^-h2^3) eieve h, la resolution 2^ denote une zone de haute frequence 
dans rimage A^j^Jk la localisation (2x,2y), En effet, les images de details etant 
obtenues par un fQtrage passe-haut de I'image A^j^J , chaque contour de A^j^J 

genfere un coefficient ondelette eieve dans une des images de detail k la resolution 
25 2^ et correspondant £l Torientation du contour. 

Ainsi, la formulation de la saillance d'un coefficient donnee dans Tfiquation 
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1 se justifie. 

S.4 Choix des points sainantsi 

Une f ois la construction des cartes de saillance achev^, nous proposons une 
m€thode afin de choisir les points les plus saillants dans Timage originale. 
5 Pour ce faire, nous construisons une arboiescence des valeurs de saillance a 

partir des -r cartes de saillance construites. De fagon analogue h la construction de 
Tarborescence des coefficients ondelettes, nous pouvons construire 2**'*^'' arbres 
de coefficients de saillance, chacun ayant pour racine un coefficient de S^r . 

Comme pour la technique Zerotree, chacun de ces coefficients correspond a une 
10 zone de taille 2x2 coefficients dans la carte iS^+i . On pent alors construire 

recursivement Tarhre dans lequel chaque noeud se voit attribuer quatre noeuds fils 
dans la carte de saillance de resolution inunddiatement sup6rieure. La figure 4 
illustre cette construction. 

Afin de localiser les points les plus saillants dans /, nous proc6dons : 

15 1 . ^ un tri par ordre decroissant des 2**^*^'^ valeurs de saillance presentes dans 

2. kla selection de la branche de saillance maximale de chacun des 2**'*^*" 
arbres ainsi tri^. 

Afin de s61ectionner cette branche, nous proposons d'effectuer un parcours 
20 de I'arbre h partir de la racine dans lequel nous s61ectionnons k chaque niveau de 
I'arbre le noeud fils ayant la valeur de saillance la plus importante (cf. figure 5). 
Nous obtenons ainsi une liste de -r valeurs de saillance : 

Branche saillante = {s^ (jc„y,),*3^ (x^,y2)X.s^^ ix_„y_,)} 
avec (x^,y^) = Arg Max{s^,^^„ (2x^.1 + u,2y^_^ + v),0 s k s 1,0 s v s l} . 

25 A partir des branches les plus saillantes de chaque arbre, le pixel de / 

choisi comme 6tant le plus repr6sentatif de la bianche est localise en (2x_^;2y_, ) . 
Dans la pratique, on parcourt uniquement un sous-ensemble (tes 2**'*^'* arbres. 
effet, pour de nombreuses applications, on recherche un nombre fixe n de points 
saillants. Dans ce cas, il convient de ne paicourir que les n arbres ayant les racines 
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les plus saillantes. 

4 Description d^aill^e d'm moins un mode particalier de 

Dans cette section, nous utilisons les 616ments techniques presents < ^ftpig la 
S section pr6c6dente dont nous fixons les pairamdtres nScessaires afin de d6crire un 
mode particulier de realisation. 

6.1 Choix de la transformation ondelette 

Comme nous Tavons mentionn^ en section 5,1, nous devons en premier 
lieu choisir une base d'ondelettes et le niveau minimal de resolution 2'' (r:^-!). 

10 Pour ce mode particulier de realisation, nous proposons d'utiliser la base de Haar 
etr=-4. 

La base de Haar est d6finie par : 



15 



\Osinon 



pour la fonction echelle, et par : 



2 

2 

0 sinon 



pour la fonction ondelette. 

^ Construct ion de Varborescence des coefficients ondelettes 

20 Dans cette etape, aucun paramfetre n*est requis. Le processus est done 

conforme ^ ce qui est d^crit dans la section 5.1. 

Construction des cartes de saiUance 

Dans cette etape, nous devons choisir les param&tres ai,(-l^k^r) 

permettant de rdgler I'importance accord^ aux coefficients de saillance suivant le 
25 niveau de resolution auquels ils appartiennent. 

Dans ce cas particulier de realisation, nous proposons d'utiliser 
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a, ^ V;tG[r,-7j 

6.4 Choir de s points^ saillantfi 
Cette etape ne requiert aucun parametre. Le processus est done conforme h, 
ce qui est decrit dans la section 5.4. 

5 §A ResuUats experimentoMx 

Les r&ultats obtenus sur des images natuielles en utilisant les paramfetres 
proposes dans ce mode particulier de r6alisation sont illustr^s sur la figure 6. 
6A Exemnle d'appUoMtion 

Parmi les applications potentielles list^s dans la section 4, nous 
10 pr6sentons dans cette section Tutilisation des points saillants pour I'indexation 
d'images fixes par le contenu. 

6.6.1 But de Vindexation damages 

L'indexation d'images par le contenu permet de retrouver, parmi une base 
de donn6es d'images, un ensemble d'images visuellement similaires k une image 
15 donnee appel6e image requite. Pour ce faire, des caract6ristiques visuelles 
(appelfes aussi descripteurs) sont extraites des images et forment la signature de 
Timage. 

Les signatures des images appartenant k la base de donn6es sont calculees 
off-line et sont stock6es dans la base de donnas. Lorsque I'utilisateur soumet par 
20 la suite une image requSte au moteur d'indexation, le moteur calcule la signature 
de rimage requSte et rapproche cette signature des signatures pr6-calcul6es de la 
base de donn^es. 

Ce rapprochement est effectue en calculant la distance entre la signature de 
rimage requSte et les signatures de la base de donn^es. Les images les plus 
25 similaires k I'images requSte sont alors celles dont la signature minimise la 
distance calculi. La figure 7 illustre ce proc€d6. 

Toute la difficult^ de Tindexation d'images consiste alors k determiner des 
descripteurs et des distances robustes. 

6.6.2 Descrivteurs has^a awr U s points saillants d'une image 

30 Nous proposons dans cette section de calculer la signature d'une image k 
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partir d'un nombre fixe de points saillants. Cette approche est inspii^ de [9]. 
Un descripteur colorimStrique et un descripteur de texture sont extraits au 
voisinage de chacun des points saillants. Le descripteur colorim6trique est 
constitu6 des moments d'ordre 0 (moyenne), 1 (variance) et 2 dans un voisinage 
.5 de taille 3x3 autour de chaque point saillant Le descripteur de texture est 
constitu6 des moments de Gabbr dans un voisinage de taille 9x9. 



entre cette signature et la signature de laf"^ image Ij dans la base de donnees est 
d6finie par : 



oil Xi et d6signent respectivement le i^'^ descripteur (par exemple i=l pour le 
descripteur colorimetrique et i=2 pour le descripteur de texture) de la image 
15 de la base et de I'image requSte R. Les poids Wj permettent de moduler 
rimportance des descripteurs les uns par rapport aux autres. 



Une fois la sigaature de Timage requdte R calculee, la distance D(RJj) 



10 
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Annexe A : Survol de la th6orie des ondelettes 

A.1 Introduction 

La th6orie des ondelettes [1][2][3] pennet d*approximer une fonction 
(courbe, surface, etc.) k difKrentes resolutions. Ainsi, cette th^orie pennet de 
5 decrire une fonction sous la forme d'une approximation grossi&ce et d'une s6rie de 
details permettant de reconstruiie parfaitement la fonction originale. 

Une telle representation multi-resolution [1] d'une fonction permet done 
d'interpreter de fagon hi6rarchique Tinformation contenue dans la fonction. Pour 
ce faire, cette information est leorganisee en un ensemble de details apparaissant k 
10 differentes resolutions. Partant d*une sequence de resolutions croissantes (rj )^ , 
les details d'une fonction k la resolution r^sont definis comme la difference 
d'information entre son approximation k la resolution et son approximation k la 
resolution ry^j . 

A.2 Notations 

t 

15 Avant de presenter plus en detail les fondements de I'analyse multi- 

resolution, nous presentons dans cette section les notations qui seront utilisees 
dans le document. 

❖ Les ensembles des nombres entiers et reels sont 
respectivement notes Zet R. 

20 *> I?(R)d6notG I'espace vectoriel des fonctions 

unidimensionnellesyfjc; mesurables et integrables. 

❖ Pour f(x)^J}(R) et g(x)^l}(R) , le produit scalaire de 
fix) et g{x) est defini par : 

(f(x).g(x)) ^ £^f(u)g(u)du . 
25 ❖ Pour f(x)^l}(R)ct g(x)E^I}(R) , la convolution 66j(x) 

et g(x) est definie par : 

/ * g(x) = fZf(u)g(x - u)du . 

❖ i^^r^^^) denote I'espace vectoriel des fonctions f(x,y) de 
deux variables mesurables et integrables. 
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❖ Pour f(x.y)E:L^(B})Qt g(x,y)E:l}(R^ le produit 
scalaiie d&f(x,y) et g(x,y) est d^fini par : 

(f(x,y),g(x.y)) -jT j^f(u.v)g(u,v)dudv. 

A.3 Propri4Ms de I'analyse multi-resolution 

5 Dans cette section, nous pr6sentons de mani&re intuitive les pTbpri6t£s 

d6sir6es de l'op6ateur pennettant I'analyse multi-resolution d'une fonction. Ces 
propri6t6s sont issues de [1]. 

Soit y4jy rop6rateur qui approxime une fonction f(x)GL^(R)h la 
resolution 2^ (J a O; (i.e.^x) est d^finie par 2^echantiUons). 
10 Les propii6t6s attendues de A^j sont les suivantes : 

1. ^yest un opdrateur lindaire. Si ^y/^jc>repr6sente 
r approximation de^fj:; k la resolution 2^' . alors A^jf(x)xie doit pas 8tre 

modifiie lorsqu'on Tapproxime de nouveau la resolution 2^ . Ce principe 
s'ecrit oA^j = A^j et montte que Toperateur A^j est un op^ateur de 
15 projection dans un espace vectoriel V^j CI?(R). Get espace vectoriel 

peut 8tre interprete comme Tensemble de toutes les approximations 
possibles k la resolution 2-'des fonctions de L^(R ) . 

2. Parmi toutes les approximations possibles de f(x) k la 
resolution 2-', A^jf(x) est la plus similaire kf(x). L'operateur ^est 

20 . done ime projection orthogonale sur V^j . 

3. L'approximation d'une fonction k la resolution 2^*^ contient 
toute rinformation necessaire pour calculer la mdme fonction k la 
resolution inferieure 2-', Cette propriete de causalite induit la relation 
suivante: 

4. L'operation d' approximation est la mSme k toutes les 
resolutions. Les espaces des fonctions approximees peuvent 6tre derives 
les uns des autres par un changement d'echelle corxespondant k la 
difference de resolution : 
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5. Lorsque Ton calcule une approximation de f[x) k la 
resolution 2^, une partie de Tinformation contenue dans^x) est perdue. 
Cependant, lorsque la resolution tend vers rinfini, la fonction approximde 
5 doit converger vers la fonctionyfj:; originale. De la m6me fa5on, lorsque la 

resolution tend vers zero, la fonction approximee contient moins 
d'information et doit converger vers zero. 

Tout espace vectoriel (V^j )^ qui satisfait I'ensemble de ces proprietes est 
appeie approximation multi-resolution de l}(R ) . 

10 A.4 Analyse multl-r4solutlon d'une fbnctlon 

unldimensionnelle 

A.4.1 Recherche d'une base de V , 

Nous avons vu dans la section A.3 que I'operateur ^approximation A^j est 
une projection orthogonale sur I'espace vectoriel V^j . Afin de caracteriser 
15 numeriquement cet operateur, nous devons trouver une base ortlionomale de V . 

V^j etant un espace vectoriel contenant les approximations de fonctions de 
I?(RJ k la resolution 2\ toute fonction f(x)GV^j p&ut gtre vue comme un 
vecteur k 2' composantes. n nous faut done trouver 1^ fonctions de base. 

Un des theor^mes principaux de la theorie des ondelettes stipule qu'il 
20 existe une fonction unique <t>(x)E:l}(R) , appeiee fonction echelle, k partir de 
laquelle on pent definir 2''fonctions de base <Pi(x)de F^^par dilatation et 
translation de 9(x) : 

^i(x) = ^(2Jx-i).i^Q,L .2^-1. 
Approximer une fonction f(x)^I}(R) a la resolution 2-' revient done k 
25 projeter de fa^on orthogonale /('jc; sur les 2-' fonctions de base ^\(x). Cette 

operation consiste k calculer le produit scalaire de f(x) avec chacune des 
2-^ fonctions de base ^{(x ) : 
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= '"^f(u).<lf(2^u^k)y^(2Ju -k). 

On peut montrer [1] que A^jf(x) peut 6tre r6duit la convolution d&f(x) 
avec le fUtre passe-bas <^(x), 6valufe au point k : 

A^jf^(f(u)*^(-2^u))(k).k^Z. 
5 Comme *!>( x) est un filtre passe-bas, A^f peut etre inteipr^t^ conune un 

filtcage passe-bas suivi d'un sous-^hantillomiage unifonne. 

A.4.2 Construction de I'analyse multi-r6solution 

Dans la pratique, les fonctions fk appioximer (signal, image, etc.) sont 
discretes. Supposons que la fonction^x; originale soit d^finie sur n = 2* (kGZ ) 
10 dchantillons. La i^solution maximale deyfx; est alors n. 

Soit 4,/ rt^proximation disa:&te 6sif(x) k la i^olution n. La propria de 
causalite (cf. section A.3) pretend que Ton peut calculer A^jfh. partir de 4,/ 
pour tout j<k. 

En effet, en calculant la projection des 2-' fonctions de base <E>/^jc ^de 
15 V^j sur V^j^i , on peut montrer que A^jf peut Stre obtenue en convoluant 
A^jufayec le mtre passe-bas correspondant k la fonction 6chelleet en sous- 
dchantillonnant le r6sultat d'un facteur 2 : 

^2jf(u) 2u)A^j^J(k),0 « < 2^ - 1 

avec h(n) = ((S>(2u),<b(u-n))yn^Z . 
20 A.4.3 La fonction detail 

Conune mentionn6 dans la propri6t6 (5) de la section A.3, 1'op^ration qui 

consiste k approximer une fonction f(x) k une r6solution 2^ k partir d'une 

approximation k la resolution 2-'*' gdn^re une perte d'information. 

Cette perte d'information est contenue dans une fonction appel^e fonction 
25 detail k la resolution 2^et not^e D^jf. n est k noter que la connaissance de 
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4j^/permetdereconstruireparfaitementlafonctionappix)xim6e A^^f . 

La fonction detail h. la resolution 2-' est obtenue en pErojetant de mani^ 
orthogonale la fonction originale/fjc; sur le complement orthogonal de V^, dans 
. Soit W^, cet espace vectoriel. 

5 Four calculer numeriquement cette projection, nous avons besoin de 

trouver une base ortbonormale de W^, c'est-Jl-dire 2''fonctions de base. Un autre 

th6or6me important de la throne des ondelettes stipule qu'k partir d'une fonction 
dchelle ^(x)^ il est possible de ddfinir 2^ fonctions de base de W^,, Ces 
fonctions de base x) sont obtenues par dilatation et translation d'une fonction 
10 ^( x) appel6e fonction ondelette : 

"^Kx) = ^(l^x - i). i = 0,L ,2-' - 1. 
De la mSme fa§on que pour la construction de Tapproximation A^,f , on 
pent monti»r que iJ^y/peut 6tre obtenu par une convolution de la fonction 
originale^A:) avec le filtte passe-haut ^^Cx) suivi d'un sous-echantillonnage d'un 
15 facteur 2^ : 

^2// - (f(u)*'9(-2^u))(k).k^Z. 

A.4.5 Extension ii I'analyse muM-rSsolutlon de fonction 
bi-dlmenslonnelles 

Dans cette section, nous pr^sentons la mani^re d'dtendre I'analyse multi- 
20 resolution par ondelettes aux fonctions de L^(R^ ) comme les images. 

Pour ce faire, on utilise les mSmes th6oremes que ceux utilises 
precedemment. Ainsi, si I'on note Tj^respace vectoriel des approximations de 
L^(R^)k la resolution 2\ on pent montrer que I'on pent trouver une base 
ortbonormale de K^^en dilatant et tiranslatant une fonction echelle 
25 ^(x,y)GL^(R^) : 

<>i(x,y) - 9(2Jx - i^Jy - J), (i.J)GZ\ 
Dans le cas particulier des approximations separables de L^(R^ on a 
^(x,y)''^(x)^(y)^^ *('jc;est une fonction echeUe de I?(R). Dans ce cas, 
I'analyse multi-resolution d'une fonction de L^(R^ Jest eifectue en ti-aitant de 
30 fafon sequentielle et separable chacune des dimensions x et y. 
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Comme dans le cas unidimensionnel, la fonction detail a la resolution 
2^ est obtenue par une projection orthogonale dej^jc,^; sur le complement de V^j 
dans V^j^i , not6 W^j . Dans le cas bidimensionnel, on pent montrer que si Ton note 
W(x)la, fonction ondelette associ6 k la fonction 6chelle 9(x), alors les trois 
5 fonctions d6finies par : 

^'(x^yJ-'i'CxJ^Cy) 

sont des fonctions ondelettes de L^(R^ ) . En dilatant et translatant ces 
trois fonctions ondelettes, on obtient une base orthonormale de FT , : 

'9f](x,y) = ^"^(l^x-kjl^y- 1 ) 

VP' (x.y) = Wra-' jc - kH^ y-l). 
10 La projection de f(x,y) sur ces trois fonctions de la base de W^j donne trois 

fonctions details : 

Dljf''f(x.y)*<^U-x)'9^(-y) 
Kf - f(x.y) * WJ(-x)^j(-y) 
D\jf'^f(x,y)*^J(-x)^j(-y) 
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revendicahons 

1. Proc6d6 de detection de points d'int^rgt dans une image num^que souice, 
ledit proc6d^ mettant en oeuvxe une transfoimation en ondelettes associant k une 
image source une image sous-6chantillonn6e, dite image €cheUe, et des 

5 coefficients ondelettes correspondant k au moins une image de details, pour au 
moins un niveau de d^»>mposition, 

un point d'int^t 6tant un point associ6 k une region de I'image pr^sentant des 
hautes frequences, 

caract6ris6 en ce qu'il comprend les Stapes suivantes : 
10 - application de ladite transformation en ondelettes k ladite image 

souice ; 

- construction d'une arborescence unique, k partir des coefficients 
ondelettes de chacune desdites images de details ; 

- selection d'au moins un point d'int^rSt, par analyse de ladite 
15 arborescence. 

2. Proc6d6 selon la revendication 1, caract6ris6 en ce qu'on determine, pour 
chaque niveau de decomposition, au moins deux images de details correspondant 
respectivement k au moins deux directions pr6determin6es par ladite 
transformation en ondelettes. 

20 3. Proc^de selon la revendication 2, caract6ris6 en ce que lesdites images de 
details conqirennent : 

- une image de details representant les hautes frequences verticales ; 

- une image de details representant les hautes frequences horizontales ; 

- une image de details representant les hautes frequences diagonales, 

25 4. Ptocede selon Tune quelconque des revendications 2 et 3, caracterise en ce 
qu'il comprend une etape de fusion des coefficients desdites images de details, de 
fagon k ne priviiegier aucune direction de ladite image source. 
5. Precede selon I'une quelconque des revendications Ik 4, caracterise en ce 
que ladite etape de construction d'une arborescence repose sur une approche de 

30 type Zerottee. 
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6. Proc6d6 selon Tune quelconque des levendications IhS, caract^6 en ce 
que chaque point de Pimage ^helle de r&olution minimum est la racine d'un 
arbre k laquelle on assode un noeud fils k chacun des coefficients ondelettes de 
chacune de la ou desdites images de d6tails localis6s au mgme emplacement, 
puis en ce qu'on associe rScursivement, k chaque noeud fils d'un niveau de 
resolution donn6, quatre noeuds fils formes par les coefficients ondelettes de 
rimage de d6tails de mdme type et du niveau de resolution precedent, et associ^e 
h la region correspondante de I'image source. 

7. Precede selon Tune quelconque des revendications 1^6, caract6ns6 en ce 
que ladite 6tape de selection met en oeuvre une etape de construction d'au moins 
une carte de saillance, affectant auxdits coefficients ondelettes une valeur de 
saillance representative de son inter§t. 

8. Procede selon la revendication 7, caract^rise en ce que I'on constant une 
carte de saillance pour chacun desdits niveaux de resolution. 

15 9. Procede selon I'une quelconque des revendications 7 et 8, caracterise en ce 
que, pour chacune desdites cartes de saillance. on fusionne pour chaque valeur de 
saillance les informations associees aux trois coefficients ondelettes correspondant 
aux trois images de detail, de fa^on h ne privildgier aucune direction dans I'image. 

10. Precede selon I'une quelconque des revendications 7^9, caracterise en ce 
20 qu'une valeur de saiUance d'un coefficient ondelette donne d'un niveau de 

resolution donne prend en compte la ou les valeurs de saillance des coefficients 
ondelettes descendant dans ladite arborescence dudit coefficient ondelette donne. 

11. Procede selon Tune quelconque des revendications 7 k 10, caracterise en ce 
qu'une valeur de saillance est une relation lineaire des coefficients ondelettes 

25 associes. 

12. Procede selon la revendication 11, caracterise en ce que la valeur de 
saiUance d'un coefficient ondelette donne est calcuiee k partir des equations 
suivantes : 
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13. Proc^ selon la revendication 12, caract€ris6 en ce que le paramfetre 
vaut -l/r pour toutes les valeurs de k. 

14. Proc€d6 selon I'une quelconque des revendications 7 k 13, caract6ris6 en ce 
5 que ladite 6tape de selection comprend une 6tape de constraction d'une 

arborescence desdites valeurs de saillance. 

15. Proc6d6 selon la revendication 14, caract&is6 en ce que ladite 6tape de 
construction d'une arborescence desdites valeurs de saillance repose sur une 
q>proclie de type Zerotiee. 

10 16. Proc6dd selon I'une quelconque des revendications 14 et 15, caract&is6 en 
ce que ladite 6tape de selection comprend les Stapes de : 

- tri par ordre ddcroissant des valeurs de saiUance de la carte de saiUance 
correspondant k la resolution minimn rn ; 

- selection de la branche pr^sentant la valeur de saillance la plus 61ev& pour 
15 chacun des arbres ainsi tri6s. 

17. Proc6d6 selon la revendication 16, caract6ris6 en ce que ladite dtape de 
selection de la branche pr6sentant la valeur de saillance la plus €\&v6& met en 
ceuvre un parcours de I'arbre correspondant H partir de sa racine. et une s^ection h. 
chaque niveau de I'arbre du noeud ffls pr6sentant la valeur de saillance la plus 

20 61ev^. 

18. Proc6d6 selon Tune quelconque des revendications 1 k 17, caract^iis6 en ce 
que ladite transformation en ondelettes met en ceuvre la base de Haar. 

19. Proc^6 selon I'une quelconque des revendications 1 k 18, caract6ris6 en ce 
que le niveau minimal de resolution est 2^. 

25 20. Procdd6 selon I'une quelconque des revendications 1 k 15, caract^ris^ en ce 
qu'il comprend une 6tape de calcul d'une signature d'image, k partir d'un nombre 
pr6ddtennin6 de points d'int6:6t de ladite image. 
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21. Proc^e selon la revendication 20, caract6is6 en ce que ladite signature est 
utUiste pour rindexation d'images par leur contenu. 

22. AppUcation du proc6d6 de detection de points d'intdrSt dans une image 
num&ique source selon Tune quelconque des revendications 1 a 21 ^ au moins un 

5 des domaines appartenant au groupe compcenant : 

- le tatouage d'images ; 

- rindexation d'images ; 

- la detection de visage(s) dans une image. 

23. Dispositif de detection de points d'int6r6t dans une image num^rique 
10 source, mettant en oeuvre une transformation en ondelettes associant h une image 

source une image sous-€chantillonn6e, dite image ^helle, et des coefficients 
ondelettes correspondant k au moins une image de details, pour au moins un 
niveau de d^omposition, 

un point d'int&et dtant un point associ6 h une region de I'image prfisentant des 
IS hautes fii^quences, 

caract^ris^ &n ce qu'il comprend : 

- des moyens d'application de ladite transformation en ondelettes k 
ladite image source ; 

- des moyens de construction d'une arborescence unique, k partir des 
20 coefficients ondelettes de chacune desdites imag^ de details ; 

- des moyens de selection d'au moins un point d'int^t, par analyse de 
ladite aiborescence. 

24. Programme d'ordinateur comprenant des instructions de code de 
progranune pour l'ex&:ution des Stapes du proc6d6 de detection de points 
d'int6r6t dans une image num^rique source selon Tune quelconque des 
revendications 1 k 22. 

25. Produit programme d'ordinateur comprenant des instructions de code de 
programme enregistr6 sur un support utilisable dans un ordinateur, 
comprenant des moyens de programmation lisible par ordinateur pour la mise en 

30 ceuvre d'une transformation en ondelettes associant k une image source une image 
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sous-6chantillonnee, dite image 6cheUe, et des coefficients ondelettes 
coirespondant k au moins une image de d^taUs, pour au moins un niveau de 
decomposition, 

un point d'int6r6t 6tant un point associd k une region de Pimage pr6sentant des 
5 hautes frequences, 

caractdrise en ce qu'il comprend : 

- des moyens de programmation lisible par ordinateur pour effectuer 
Tapplication de ladite transformation en ondelettes h ladite image 
source ; 

10 " des moyens de programmation lisible par ordinateur pour effectuer la 

construction d'une arborescence unique, h partir des coefficients 
ondelettes de chacune desdites images de details ; 

- des moyens de programmation lisible par ordinateur pour effectuer la 
selection d'au moins un point d'int6r8t, par analyse de ladite 

15 arborescence. 

26. Support de donnSes num6riques utilisable par im ordinateur, caract6rise en 
ce qu'il comprend des instructions de code de programme d*un programme 
d'ordinateur selon Tune quelconque des revendications 24 et 25. 
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